x <- c(6)
y <- c(2)
alpha <- 1.47E-2
beta <- 1.63E-3
gamma <- 1.2E-2
delta <- 3E-4

for (i in 1:4000) {
  lastx <- x[length(x)]
  lasty <- y[length(y)]
  dx = lastx*(alpha - beta*lasty)
  dy = -lasty*(gamma - delta*lastx)
  x <- c(x, lastx + dx)
  y <- c(y, lasty + dy)
}

#plot(x, type="l", col="green",ylim=c(0,max(x)))
#lines(y, type="l", col="black")

plot(x,y, type="l", col="gray", xlab="pavi\u0101ni",ylab="gepardi")


